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Abstract. We give the solution of certain parabolic evolution problems (time-depending 
perturbations of the heat equation for the harmonic oscillator ) as explicit integrals on 
a space of continuous functions, called the Wiener space. The methods are based on 
the Mchlcr formula giving the solution of the unperturbed problem and on the use of 
discretization to split the difficulties. 



The aim of this article is to give an explicit expression for the solutions of the problem 



as a Feynman Kac type integral on the Wiener space Cw, which is the space of continu- 
ous functions on [0, 1], equipped with a special measure mw called Wiener measure (see 
Subsection 12. ip . 

Recall that the theory of semigroups and their perturbations gives the existence and 
uniqueness of the solution of ([I]) in C°(R + , L 2 (R)), if, for example, c is continuous, 
bounded and satisfies the following Holder condition 

(2) 3L > 0,3a G [0, 1] : Vs,i G [0,oo), Va; G R, \c{t,x) - c{s,x)\ < L\t - s\ a . 

The problem of giving such explicit solutions has already been studied under more re- 
strictive conditions than in the present case. In [3], [6], the case when the potential 
V(t, x) = x 2 + c(t, x) is replaced by a general, unbounded but time-independant potential 
V(x) is studied and in [5], the author deals with the case when V(t,x) is bounded. 
The present work is devoted to the study of the mixed case ([I]) . The main result of this 
paper is the following 

Theorem 1. Let Vq G L 2 (R). Let c be continuous and bounded on [0, oo) x R. Suppose 
that c belongs to L 2 (]0,T[xR) for any T > and satisfies the Holder condition |D). The 



1. Introduction. 



(1) 




(t,x) + {x 2 + c(t,x))v(t,x) = on ]0, oo[xR 
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function v defined on ]0, oo) x R by 



v(t, x) 

(3) 



J v (V2t w(l) + x) exp (^-t J (x + V2tw(s)^ ds^j 

exp ^— t J c (t(l — s), y/2t w(s) + x^j ds^j dmw(w) 



is 



the solution of (QJ) in C°(R + , L 2 (R)). 



This result differs from the usual Feynman-Kac formula, for the function w G Cw is 
never evaluated in t. No well known change of measure ([6], PQ, [9]) on Cw allows to get 
one formula from the other one. 

One interest would be to get regularity results (or to study the dependency on a param- 
eter) more easily than using the theory of semi-groups, by straightforward derivations. 
Problem (pQ) will be considered as a perturbation of the heat equation related to the har- 
monic oscillator H. Section[2]contains the notions needed about Wiener integrals, followed 
by some facts about the heat kernel for the harmonic oscillator, in particular Mehler's 
formula. The proof of a preliminary version of Theorem [H under restrictive regularity 
conditions, is completed in Sections [3] and 0] by constructing a sequence of functions ex- 
pressed as integrals on Cw and converging to the solution of ([!]). The regularity conditions 
are weakened in Section [6], which leads to Theorem [TJ In section a slight modification 
of section [3] gives an alternate demonstration of Feynman-Kac formula, which does not 
rely on Ito calculus. Eventually we give some explicit expressions of Wiener integrals, 
which can be deduced from Theorem [T] ( annexe . 

2. Preliminaries 

2.1. Wiener integrals. The construction of the Wiener measure is detailed, for example, 
in [6] and [ID] . The (classical) Wiener space, denoted by Cw, is the set of all real- valued 
continuous functions w on [0, 1], with u>(0) = 0. It can be equipped with a probability 
measure mw called Wiener measure, which is defined on a a-field 9JT* containing all sets 
of the type 

J = {weC w ■ (w(ti),...,w(t n )) G U}, 

where to = < t\ < ■ ■ ■ < t n < 1 and % is a Borel set of W 1 . For this kind of set the 
measure is given by 

m w (J) = / /ti,...,tn 60 d£i 

Ju 

where ft lr ..,t n is the normal density 




//. o (airlift-** 

V i=i 

(with £ = 0). 

The integral of an integrable, real-valued function F defined on Cw will be denoted by 
or f Cw F(w) dmwiw). When F(w) depends on the values of w at finitely many 
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fixed points < h < t 2 < . . . t n of [0, 1], E W [F] can be expressed as follows. Let (p be 
a measurable real-valued function on R n and let 1Z 1 denote the Borel sets on R. The 
mapping 

F : (C W ,M*) ->• (R,^ 1 ) 

u> 1 ^ ip(w(tx), ... ,w(t n )) 

is measurable and 



(4) / F(w)dm w (w)= <p(£i,'",€ n )ft u ...,t n (£i,'-,€ n ) dfi...dfn, 

in the sense that the existence of one side implies the existence of the other and the 
equality. When F depends on the value of w at infinitely many values of t, it is useful 
to recall that, if the topology on Cw is defined by the uniform norm, any open set is 
9JT*-measurable and any real-valued, continuous function is OJT-measurable. 
One important propertie is that the so-called coordinate process, defined on Cw by 

Vt G [0, 1], Vw e C w , X t (w) = w(t) 

is a Brownian motion. A consequence is Fernique's Theorem, which states that there 

e*,sts a positive „ such that the integra. / e*»' dmw(w ) converges, || u || be,n g the 

uniform norm of w (|6]). 

2.2. Heat equation for the harmonic oscillator. Classically the problem 

° V (t, x) + x 2 v (t, x) = (t, x) e]0, oo[xR, 




(5) 

(where t>o G L 2 (M)) has a unique solution v belonging to C°([0, oof, L 2 (M)) ([8J). More 
precisely, there exists a semigroup (U t )t>o of L 2 (M)-contractions such that, for all t > 0, 
f (t, •) = C/tDo- In the case of the harmonic oscillator this semigroup is explicitely given by 
Mehler's formula ([!]) 

One can check directly that v : (t, x) 1— > (UtVo)(x) is infinitely derivable on ]0, oo[xM. The 

change of variable z — - = y gives another expression, 

ch(2t) 

(7) v(t,x) = (U t v )(x) = / q(t,x,y) v (y+ X ) dy 



ch(2t) 



where 



Both formulae give solutions of (E]) for initial conditions which do not necessarily belong 
to L 2 . For example, if vq is continuous and vanishes at infinity, U t Vo is C°° on ]0, 00 [xK. 
and continuous on [0, +00) x R. The property 

Vt, r > , U t U T v = U t+T v 
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remains valid as well. 



3. Approximation sequence 

This Section is the first step in the proof of Theorem [TJ It is devoted to the construction 
of a sequence (v„ ) nS N of functions designed to approximate the solution of flTJ). For a 
given positive integer n and a fixed positive t, n factors e^ H coming from (jSJ) alternate 
with n factors en c ( tk > x ) (as in Trotter's formula), where the tk s are the bounds of a suitable 
discretization of [0, t]. The expression thus found can be written as an integral on the 
Wiener space which converges to the function v in Theorem HJ The exponent (t) stresses 
the fact that the sequence depends on the upper bound of the time interval. 

Suppose that the function c = c(t, x) is continuous and bounded on [0, oo) x R and that 
v o — v o{ x ) is continuous on K and vanishes at infinity. Let us split the interval [0,t] into 
2n subintervals bounded by the = ||, k — 0, . . . , 2n and assign to satisfy, on each 
subinterval, an incomplete version of 

3v d 2 v 

— (t, x) - Q^{t, x) + x 2 v(t, x) + c(t, x)v(t, x) = 0. 

The sequence Vn^ is constructed by induction in the following way 

• on the even interval [T 2 fc, r 2 fc + i], v$ is a solution of the heat equation for the 
harmonic oscillator 

dv , N _ / d 2 v , N 2 



and satisfies the initial condition lim v(r,x) = v^(T 2 k,x), where vfflfak,') was 

T — > T 2k 

built in the preceding step ; 
• on the odd interval [r 2 k+i, T2&+2], Vn* satisfies the ordinary differential equation 

dv 

(9) — (r, x) + 2c(r 2 fc+ 2 , x)v(r, x) = 0, 

or 

with the initial condition lim v(r,x) = v^\r 2 k+i,x), where Vn\r 2 k+i,-) was 
built in the preceding step. 
Both equations have of course an explicit solution 

v ( n\r,x) = {U 2i r-r 2k )V^\T 2k ,-)) (x) on [r 2k , T 2k+1 ] , 

v®(t,x) = e- 2c ^^ T - T ^v^(T 2k+ ux) on [r 2fc+1 , r 2k+2 \. 

The equation (Q is a constant coefficient linear differential equation since c depends on 
the fixed time T 2 k+ 2 - The initial conditions ensure the continuity with respect to r. The 
factor 2 in OH]) and ()9]) compensates the fact that each equation is solved on half of the 
interval (see Lemma |9] below and [5], where this discretization was introduced, to treat 
the usual heat equation.) 



A REMARK CONCERNING FEYNMAN KAC FORMULAS FOR THE PERTURBED HARMONIC OSCILLATORS 
.(*)/ 



We first write Vn (t,x) as iterated integrals. 
Proposition 2. For all real x we have 



vS\t,x)= [ (2nsh(2t/n)r n/2 v (o- n 

/ 1 ch(2t/ra) ^ / a. 
exp -o^ToTT^Z. 



X 

+ 



ch(2t/ra) 

n-j 



exp 



2sh(2t/n)^V ch(2t/n) 
1 sh(2i/ra) ^ / x x " 

2ch(2t/n) ^ V n_j + ch(2t/n)"-J 



n 

exp 



V n ' ch(2t/n) n -^ 
wjt/i £/ie convention a — 0. 

Proof. By induction on fc we get that, for all r G [r 2 fe, T 2 fe + i], 
v®(t,x)= q(2(T-T 2k ),x,y k+1 ) 



7»% 



II 9 I 5Z c h(2t/n)'- 1 -i + ch(4(r - r 2fc ))ch(2t/n) fc - 

(10) fc / * / 

JJexp ( --c ( r 2j , J] ch (2t/ n )n^' + ch(4(r - r 2fc ))ch(2t/n) fe - ; 

/fe+1 \ 

V ° C h(2t/n)^ + ch(4(r - r^ch^/n)* J d fi ■ ■ ■ 
For r = r 2n = t and A; + 1 = n, this gives 



n V 1=1 



ch(2t/n)^ 1 ch(2t/n) n 



11^ | S C h(2tjny-^ + ch(2t/n)^" yj 



j=i \ i=j+i 



n 1.1 n 



I! ex P ; E eh(2t/n)'-W + M2t X /n)^ • • ^ ' ' ' ^ 
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with the convention that a sum is equal to zero if its lower index is strictly superior to its 
upper index (this is useful for j = n). Let 



n 



(H) ak= £ d#p^ , 1 < k <n 

l =n +l-k K 1 ' 

(recall a — 0). Reciprocally, 

and this change of variables y — > a has Jacobian equal to 1. Hence 

W?( '' X) = W*" + «^) 

ft q (n' + (ch(2t/n))-J' Gn+1 - j ~ ch(2t/n)' ) 

A ex P ("^ ^ + (ch(2 J n)) ^ ) ) ^ • • ■ ^ 

Replacing q by its expression gives the Proposition. □ 

The second step is the transformation of this integral on M™ into an integral on the 
Wiener space. Since CV does not depend on n we then shall be able to let n converge to 
infinity. These computations lead to the following Proposition. 

Proposition 3. Let v be continuous and bounded on ffiL Suppose c is continuous and 
has a lower bound on [0, oo) x R. Moreover let c have an x-derivative c x at each point 
of ]0, oo) x K. and suppose c x has an extension to [0, oo) x R which is continuous and 
bounded. 

When n goes to infinity the numerical sequence Vn\t, x) converges to a limit v(t, x) defined 
by 

v(t,x) = j v (V2t w(l) +x) exp^-t^ (x + V2tw(sfj ds^j 

exp ^— t J c (t(l — s), V2tw(s) + x^j ds^j dmw(w). 

Proof. To bring out the normal density one uses the change of variables 
(12) £j = <7> sh(2t/n))- 1 / 2 , 1 < j < n 
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and the convention £ = °o = 0. It follows 

v®(t,x) = [ (2n/n)- n / 2 v ((nsh(2t/n)y/%+ ' 



ch(2t/n) n 



exp - 



X 

n) 



exp 



We introduce the normal density for equidistant points < U — i/n < 1 , 

f tl ,...,t n (tu---^n) = (^/nr n/2 e W 
to make the transition to Wiener space more natural. We first obtain 



v { :\t,x) = I v ((n sh(2t/n))V^ n + x ) ./;, lA r : . 



ch(2t/n) 7 



( n ^ ™ / 

- - e«-i) 2 - ^(^ ch (^/^) E 
i=l 3=1 \ 



t \ 2^ 



ch(2t/ra) 



exp 



n— 1 / \ 2 N 

X 



2 ch(2t/n) ^ V ch(2t/n)J 



3=0 

and then, by formula (H 



= / F n («;)G n («;)ir n («;) ^ ( (n sh(2t/n)) 1 / 2 «;(l) + X ) dm w (w), 
Jc w V cn^r/n; 
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with 



F n (w) = exp ^ [(1 - ch(2t/n)) W (t i ) 2 + (1 - ch(2t/n)- 1 ) W (t i _ 1 ) 2 ] 



i=l 



*„(«) = exp (-1 gc („ ^a/nM^O/n) + 

The Proposition is now a consequence of the dominated convergence Theorem. The 
convergences and estimations concerning the four factors are treated in the Lemma just 
below 

Lemma 4. Denote by m a lower bound of c. For all w G Cw, 

• < F n {w) < 1 and lim F n {w) = 1, 

• < G n (w) < 1 and lim G n (w) = exp ^— t J (^x + \Flt w(s)j ds*j , 

• < H n (w) < e~ mt and lim HJw) = exp i-t [ c (t(l - s), \Plt w(s) + x) ds), 

n->oo JO ^ ' 



v ( (n sh(2t/n)) 1/2 w(l) + 



;r 



lim w ( (n sh(2t/n)) 1/2 w(l) + 



ch(2t/n) n 



ch(2t/n) 



< ||fo||oo OWcZ 

= ^(v 7 ^ w(l) 



Proof of Lemma [7] 

We shall use more than once the continuity of w G Cw to compute limits of Riemann 
sums like lim^oo X]j=o w (^j)- The positivity and the estimates proposed are obvious, 
except for F n which we shall treat first. 
Let us develop the argument of the exponential 

n 

A n = | X) K 1 - ch(2t/n))w(U) 2 + (1 - ch(2t/n)- 1 ) W (t l _ 1 ) 2 ] . 
i=i 



Since t n = 1 and w(t ) = w(0) = 0, 



n /. ,,„,,, ,„ s9 n (ch(2t/n) — l) 2 v-^ / x9 
A n = -(1 - di{2t/n))w{lf - — - — ~n7oT7T\ — ^ ' 

j=i 



2 ch(2£/n) 



and then F n < 1. 

Now we claim that A n converges to 0. An asymptotic expansion is enough for the term 
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containing w(l). The second term can be written as 

n 2 (ch(2t/n) - l) 2 1 " 
~2 ch(2t/ra) X n , 



In this product the second factor converges to J* w 2 ds and the first one to 0, using an 
asymptotic expansion. 

To compute the limit of G n we split the argument of the exponential into three terms : 



B„ = 



71 — 1 

X 



S -^M-y(nsh(2t/n)) 1 / 2 w(t J ) , 
ch(2t/n) V 1 JJ KjJ ch(2t/ny 



1 sh(2t/n) x ^ x- 
n ~ ~2ch(2t/n) ch(2t/n) 2 ^' ' 

A direct computation gives 

n 2 sh(2t/n) 2 1 £=i , x2 (2t) 2 f 1 2/ , , 

A n = \ 7 ; ; -> w(tj) — > / w s ) ds . 

2 ch(2t/n) n 1 j; 2 ] U 

The third term is a geometric sum. An asymptotic expansion leads to 

Cri y tx . 

The second term is decomposed as 

l «' 1=0 V « 7 7=0 



rt— 1 



Since 



one has 



Moreover 



n 

j=o j= " 

v ' V 

D n E n 



n 3 'hh(2t/n) 3 / 2 ,„ 

-w =(21) (1+ ° (1/n)) ' 



L>„ — >■ -x(2t) 3/2 f w(s) ds. 

Jo 



1 1 

ch(2t/n)i ~ ~ ch{2t/n) n 



< 1 - , ^ . . < 1 - T . . r 2t 2 T/T 1 
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^1*' ch(Wn) (1 -* (8/n) )} ng K ^' 



0. 



To sum up, (ch(2t/n)) J can be replaced by 1 in all the w(tj)(ch(2t/n)) j . We conclude 
that 



+ B n + C n 



-2t 2 I w 2 (s) ds -x(2t) 3/2 [ w(s) ds - tx 2 
Jo Jo 

-t / (x + {2tf/ 2 w{s)) 2 ds. 
Jo 



Let us turn to H n . For all j and all w we can write 

_ c| (^, ( „ sh(2t/n)) v 2m o /n) + _|_ 



< —to, 



hence the estimate H n (w) < e mt . Denote by | \c x \ |oo the uniform norm of c x on e [0, oo) x 
R. The mean value theorem implies that 



c( t -i t ,(« S h(2 t /n))V Mj /«) + ch(2t/n)J 



c (t - -t,V2t w(j/n) + x 



'n sh(2t/n)) 1/2 w(j'/n) - v 7 ^ w(j/n) 



+ C x 



^ Qc oo P oo 



(n sh(2t/n)) 1 / 2 - V2t 



Cx oo 



X 



ch(2t/n) n 



ch(2t/ny 
- 1 



x 



in which the last term is independant of j. Therefore 

n— 1 . ,-, , n— 1 



71 



£c(f - J -t, (n shA^Mj/n) + — 1-) ^ v^t «;07n) + s) 



j=0 

< illCxIloodklloo 



n 



sh(2t/n)) 1/2 - V2t +\x\ 



3=0 



I 

n 



- 1 



). 



ch(2t/n) n 

which converges to when n goes to infinity. We deduce that H n {w) and 



n-l 



cxp 



3=0 



-VcU-H v 7 ^ w(j'/n) + x 
n { V n 



converge to the same limit. Since the function u ^ c (t — ut, y/2t w(u) + x) is continuous 
on [0,1], this limit is 



exp ^— t J c (t(l — s), \p2± w(s) + x^j ds^j 
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Let us treat the last point of the Lemma. Recall that Vo is continuous and bounded. 
Its argument goes to {2t) l / 2 w{l) + x because 

(It 2 1 \ 
ch(2t/n) n = exp — + o(— ) ->• 1, 
\ n n 2 J 

which completes the proof. □ 



4. Preliminary version of Theorem [TJ 

We still need to show that the function v constructed above as the limit, at time t, of 
the sequence (v$) ne n is a solution of <^j. The demonstration requires stronger regularity 
conditions on vq and c than the ones used to compute the limit. Here is the ( weaker ) 
version of Theorem [T] which will be proved in this section. 

Theorem 5. Let vq be a C 4 function over R 7 which has bounded derivatives of order up 
to 4. Suppose Vq(x) converges to when x goes to infinity. 
Let c be a function which 

• is continuous and bounded on [0, oo) x R, 

• is C 1 on ]0, oo) x R, 

• has bounded space derivatives, up to order 4, these derivatives being continuous 
and bounded on ]0, oo) x R. 

The function v defined on ]0, oo) x R by 
v(t, x) 



j v (V2t w(l) +x) exp y (x + V2tw(s)^ ds 

exp ^—t J c (t(l — s), V2t w(s) + ds^j dmw{w) 

is a solution of (T7]j. 

We need preliminary results. The first Lemma shows that the sequence (vn ' > ) and one of 
its derivatives converge on a dyadic subset of [0, t}. The second one gives uniform estimates 
concerning some of the derivatives of v„ . The last result proves that a subsequence of 
(v$) converges uniformly, and on [0,t] itself. 

Lemma 6. Let t g]0, oo[ ; let T> be the set 

k 



V 



[re[0,t] : 3n eN,3A;oe{0,...,2" -l},r=^t} 



For all r GD and all iGl, 

(t), \ /x , d 2 Vo} . . d 2 v . , 

lim v^i{T,x) = v{T,x) and hm - (t,x) = -7t^(t, x) 

n-^oo n-s>oo OX OX 

where v is the function defined in Theorem^ 
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The rest of this Section is devoted to the proof of the Lemmas and of Theorem [5j 
Proof of Lemma 

It is easier to express Vn\r, x) when r = kt/n is a bound of the subdivision. Therefore we 

k 

are led to consider nested subdivisions. A dyadic point r = —t, which is already a bound 
of the subdivision with 2-2 n intervals, is a bound of all following subdivisions. The set T> is 
the unions of all such points. Formula (fTUl) gives Vn\r, x) for r G [-z~t, — t). Adapting 

—,11 / ' 

k 2 p k 

the proof of Proposition ([3]) we deduce the expression of v^„ +p (t, x) when r = —t = ^ n+p t 

and obtain 
(13) 

x 



^ x ) = J v ° sh(2t/2^) W (l) + ch(2t/2w+p)2Pfc 



2Pfc-l , / , \ 



exp 



lsh(2t/2"+f) V^V ,— — / 



2Pk ) ch(2t/2 n +P)' 




iTPk. . _ . x2 2 p A; (ch(2t/2 n+ P 

exp — (1 - ch(2t/2 n+p )w(l) 2 - v v 1 



2 v v ' y w 2 ch(2t/2 



n+p 




(I 2 Vt^ 2 p k — I I x \ 
> c( : — t, j2Pk sh(2t/2 n +P)w(—) + — — ; — — ) dm w (w). 
2«+p ^ 2 n+ P v v ' ' y 2Pk' ch(2t/2 n+ P) lJ I v ' 

The integrated terms are similar to the F n , G n , H n treated in Lemma HI We get the limit 
of v^} +p (t,x) by letting p go to infinity in the integral and do not need the additional 
hypotheses. The only difference is that 

lim 2 p k sh(2t/2 n+p ) = r. 

To find the limit of 9^i4«(t, x), one has to derivate ({TBI twice with respect to x. The 
derivatives of the first exponential term contain the expression 

sh(2f/2»+») / / ( \ , x \ 1 



ch(2t/2"+P) ^ V V 2 ^/ ch(2t/2 n +f)7 



ch(2t/2 n +P) r 



as well as its square and its derivative (with respect to x). We estimate M by K(\ |u>||oo + 
\x\) with a constant depending only on t. The other exponential terms and their deriva- 
tives are bounded with respect to w. Therefore we can apply the dominated convergence 
Theorem, since f c \\w\\ 2 dm w {w) is bounded (by Fernique's Theorem, see Section [2]). 
As for the limits themselves, the same techniques can be applied as in the proof of Lemma 
HI Note that the estimate of d^c is needed to treat terms containing d 2 c for we use the 
mean value theorem. □ 
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Lemma 7. There exists C > 0, depending on t and x, such that 

~kt (k + l)t 



Vj G {0, . . . , 4}, Vn G N*, VJfe G {0, . . . , 7i - 1}, Vr G 



n n 



d j v$ (r, x) 



< C. 



Proof of Lemma [7| 

We first treat the interval [f^, Formula ([TDD shows that 



V n (T,X) 



2n ' 2n 

(27rsh(4(r - kt/n)))- 1/2 (2irsh(2t/n)y k/2 



exp 



exp 



exp 



1 sh(4(Y - kt/n)) 
2ch(4(r - kt/nj)' 



x exp 



1 ch(4(r - kt/n)) 2 
2sh(4(r- kt/n)) ai 



lch(2t/n) 
2sh(2t/n) ^ 

\ ' ' 1=1 



+2-./ 



ch(2t/n) 



lsh(2£/n) 
"2ch(2t/n) 



7=1 ^ 



X 



ch(4(r - kt/n))ch(2t/n) k -i 



exp 



c ( -t, o- k+1 -j + 



x 



ch(4(r - kt/n))ch.{2t/n) k -i 



x 



da 1 . . . da k+1 . 



ch(4(r - kt/n))ch(2t/n) k / 

The space-derivatives of order at most 4 of all terms but one are bounded by constants 
depending on t, \\vq , || and ||9|c|| (j < 4), but not on k, n, r or x. The derivative of 



1 sh 
2ch ( 



n 



exp ( ---f 1 ! KV - — )).r- ) exp [ ( + 



ch(4(r-f))ch(f)^-i 
is less easy to treat. It is the product of the exponential term and of 

1 



sh . . kt . , sh . 2t . v-^ / 



77 



J = l 



ch(4(r - f ))ch(f y ch(4(r - f ))ch(f ' 



We then have to estimate expressions such as 



(^2 a s b s c s y exp a s ^J 
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where the a s are quotients sh/ch and the c s are the (ch(4(r — — )) -1 ch(— ) _s . Applying 
Cauchy-Schwartz inequality to the sum outside of the exponential term shows it is smaller 
than \fYl< a sb 2 s 2~2 a sC 2 s - The first factor is absorbed by the exponential and the second one 
is a geometric sum, which is bounded. 

It is eventually possible to estimate the derivatives of order < j < 4 by 

C / (2vrsh(4(r - kt/n)))-^(2n S h(2t/n))-^ exp ~_ j^fo 

lch(2t/n) A/ (Tk+i-j \\j A 

exp i -2 iu^M g r +2 - j ~ J J 1 • • • k+i ' 

where the constant C depends on t, \\vq \\oo et H^'cHoo (0 < j < 4) only. The integral is 
equal to 

-1/2 



chf4(r-— )J ch{2t/n)~ k / 2 < 1, 



which establishes our claim for the interval ]. 
On the following interval, [ 2 ^-t, 2|±?], 

„J)(t,s) = exp (-2(r - ^t)c(^i) lZ )) «^^*.*)- 

It is the product of two functions having bounded space-derivatives (of order at most 4). 
This proves the estimations on [0, £]. □ 

Lemma 8. Let {u n ) n& n be a sequence of functions continuous on [0, T], piecewise C 1 on 
[0, T] and satisfying 

3C > : Vn € N, Vr G [0,T] swc/i t/iat <(r) exists, \u' n {r)\ < C. 

Suppose D is a dense subset of [0,T] and u, a continuous function on [0, T], such that 

Vr G -D, limw n (r) = w(r). 

n 

Then the sequence (u n ) has a subsequence which converges to u uniformly on [0,T]. 

Proof. The bounds on the derivatives show that (w n ) ng N is bounded and equicontinuous 
on [0, T]. Ascoli's Theorem yields the existence of a uniformly converging subsequence. 
Its limit u is continuous on [0,T] and equal to u on the dense subset D, which justifies 
the Lemma. □ 

Proof of Theorem 

Let x be fixed. We have proved that the (numerical) sequences (v® (t, x)) n£ N and {{d^v® (t, x)) n£ N 
converge (respectively) to t>(r, x) and d 2 v{r, x) for all r belonging to the dense subset V 
defined in Lemma |6j Now to prove the uniform convergence on [0, t] we apply Lemma |HJ 

The continuity of ?4n(-,x) is a consequence of its definition. The continuity of c^-i^n (•, x) 
can be proved by studying the expressions appearing in the proof of Lemma [7] and so is 
the derivability with respect to r. 

To establish the bounds on the r-derivatives let us recall that 



A REMARK CONCERNING FEYNMAN KAC FORMULAS FOR THE PERTURBED HARMONIC OSCILLATOR5 

(1) On the even intervals [2H/2.2 n , (2k + l)t/2.2 n ], 



dv$(r,x) _ n d 2 v$(T, x) 2 (t) 



dr dx 2 
(2) on the odd intervals [(2k + l)i/2.2 n , (2k + 2)t/2.2 

dv^l (r, x) 



- 2x v\l (r, x), 



Or 



-2c((2k + 2)t/2.2 n , x)v$ (r, x). 



The bounds on the time derivatives come from the estimates on the space derivatives : to 
show that drV^l is bounded one needs the x-derivatives up to order 2, to treat <9 T (<9 x ) 2 t>2^ (r, x) 
one needs the x- derivatives up to order 4. The bounds concerning the space-derivatives 
have been established in Lemma CD It follows that a subsequence of (v^l) (resp. of (d T V2n) 
) converges uniformly on [0,t]. Let us denote its indexes by (p(n). 

Next we write more concisely the system of equations defining Vn (according to the 
interval). On [0,t] \ {kt/2 n }, < k < 2 n , v$ satisfies 



dv n ) (r, x) 



with 



Wn(r) 



Pn(r) 



d 2 v%'(T,x) 



dx 2 



x 2 v®(t,x) - 2(1 - /3 n (r))c n (r,x)^(r,x), 



on 



on 



2k 2k +1 

— t, t 

2n 2n 



2fc+ 1 2fc + 2 ' 
■t, 1 



2n 



2n 



and c n (r, x) 



'2k+2 



2n ' 



t, x) on 



2A; 2A; + 2 

— t, 1 

2n 2n 



This equation still holds for the subsequence indexed by <p(n). The uniform convergence 
allows us to integrate on any sub interval [0, s] of [0,t] : 



V l\n)( S i X ) -<(n)(M 



(t) 



9 2 4l } (r,x) 



P<p(n) (T ) 



dx' 



x v 



n)( T ^) 



/3 v (n)(r))c v(n) (r, x)t»2 n) (r, x) dr. 



if(n) 



Now let n tend to infinity. The following result ([5]) shows what become of (3 n and of the 
factor 2 : 

Lemma 9. Let (3 n be the function defined above. Let (ip n ) be a sequence of functions 
belonging to G L 1 ([0,t]) and suppose it converges uniformly on [0,t] to a limit ip. Then, 
for allO <t < a < t, 



lim / (3 n ip n ds = lim 



f3 n )ip n ds 



ip ds. 
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Eventually, for all s < t, we obtain 

v (s, x) — Vq(x) — [ ( ~g 2^ —x 2 vt,x) ) — <-•( r. .<•)<■( r. ./•) o'V . 



o 



As does not depend on t (the intermediates v$ depend on t but not the limit), 

the function v is a solution of (P). □ 

5. Another proof of Feynman Kac formula 

Before completing the proof of Theorem [1] we shall see that a small modification of the 
method developped in Section^] yields the following expression for v. For sufficiently small 

t, 

v(t,x)= / v (w(2t) + x) exp I — - / (x + w(s)) 2 ds 
Jc w \ * Jo 

( 1 f 2t \ 
exp I — - J c (t — s/2, w(s) + x) ds J dmw{w) . 

With u(t,x) = v(t/2,x), u satisfies 

du 1 d 2 u 1 

dt^ X) ~ 2 dx^ (t ' ^ + 2^ + C(t/2 ' X ^ U ^ X) = °' 
which explains the differences with the usual expression. 

The first point is that the demonstration does not use the Ito integral at all. What 
is, perhaps, more significant is that both expressions of v are not linked by a "classical" 
change of variable on Wiener space. 

Proof. It essentially follows the same steps as in Section |3j Starting from Proposition [21 
we consider the time sequence (tfc)o<fc<n with = k— instead of = k/n. This yields 

^(t,x)=f LJo^ V^o (w(2t) 



nsh(2t/n) ) \ y ' ch(2t/n 



exp 



1 ch(2t/n) ^ 
2sh(2t/n)^ 



w(2jt/n) 



w(2(j - l)t/n) 
ch(2t/n) 



2 N 



n 

n 



ex P ( 77 £ k(2i*/") - «W - l)t/n)f 



exp 



n— 1 / x 2 N 

X 



1 sh(2t/ra) s-^f fn . . s 



2ch(2t/n)^V. ch(2t/n)J 



(i n_1 fin — j)t x \ \ 

-n]^ C {-^ M2jt/n) + d^w)) dmw{W) - 
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Most of the terms are bounded and converge as in Section [3] or even more easily. It just 
remains to treat 



lch(2*/n) 
2sh(2t/n) ^ 



E 



w(2jt/n) 



w(2(j - l)f/n) 
ch(2t/n) 



+ -J2MW/n)-w(2(j-l)t/n)f 

3=1 



This expression splits into A n + B n where 



n lch(2t/n)\ , x9 \X . . x9 ch (2tn - 1 

-4n = I w / n w 2i ~ > w(2jt n) 2 \ / //, , 

V 4t 2sh(2t/n) ' ^ ' o.oU'Hl-n^CH 



3=1 



2ch(2t/n)sh(2i/n) 



and 



1 \ n_1 



Z?„ = ( . i > ic(2it/n) (ir(2itln) - w(2(j - l)t/n)) . 

The first term is negative and converges to 0. The second one can be estimated as follows 

\BJ < 



n 



1 



2t sh(2t/n) 



n 
2~t 



n 



1 



2t 



n-l 



- (w(2jt/n) - w(2(j - l)t/n)f 



3=1 



< My/2t\\w\ 



\ 



2t 



n 



n—1 



J2(™(W/n)-w(2(j-l)t/n)f 



3=1 



As J^jLi (w(2jt/n) — w(2(j — l)t/n)) is the quadratic variation of a Brownian motion, 

the subsequence for n = 2 P converges to \/2~i and it is smaller than 4rz||w|| 2 . To sum up, 
\B n \ < 4tM||w|| 2 and converges to 0. This estimation and Fernique's theorem allow us to 
use Lebesgue dominated convergence theorem, provided t is small enough. □ 

6. Proof of Theorem Q] 

To get the optimal form of the Theorem, it remains to prove that formula fl3]) gives 
a solution of Problem ([1]) even if Vq and c satisfy much weaker assumptions. This will 
be done by approximating general vq and c by regular functions and showing that the 
solution of the approximating problem converges to that of the real problem. 

Proposition 10. For vq G L 2 (M) and c measurable and inferiorly bounded on ]0, oo) x K. 
we define, following formula (T3)) ; 



S(v ,c)(t,x) 



J v {\/2i w(l) + x) exp^-t^ (x + V2tw(s)y ds^j 

exp f—t J c (t(l — s), \p2t w(s) + x^j ds^j dm w {w) 
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For all < a < j3 < oo, S(vo,c) belongs to L 2 ([a,(3] x M). Moreover, 
(14) I \S(v ,c)(t,x)\ 2 dtdx < \\v \\ 2 [ e -*™V) dt . 

Proof . 

Let us consider 

If \v (V2t w(l) + x)| 2 exp i-lt ! (x + V2tw(s)) 2 ds] 

a JrJc w V JO ^ ' J 

The first exponential factor is smaller than 1 and the second one, than exp(— It inf (c)). 
By Fubini's Theorem, 

/ < / e -2* inf ( c ) f (t \v (V2t w(l) + x)\ 2 dx] dm w {w)dt 

J a Jc w V-'K / 

< [\- 2ti ^ [ ( [ MOI 2 ^ dm w (w)dt, 



Ja JC W \JR / 

with £ = y/2t w(l) + x. Eventually, 

I< \\vo\\ 2 [\- 2ti ^ dt. 

This shows that the integral 

\v (V2tw(l) + x )|V 2t ^ 1 (-+v / 2*^)) 2 ds e -2tf \(t(i-s),Vn „(.)+*) ds dmw{w) 

Cw 

converges for almost all (t, x) G [a, 0\ x R. By Holder's inequality it follows that 
S(vo,c)(t,x) is defined for the same (t,x) and that S(vq,c) satisfies the inequality (fT4j) . 
□ 

Proposition 11. Let (v^) G L 2 (R) N converge to v in L 2 (R) and (c (n) ) G L 2 (M+ x M) N 
converge to c G L 2 (IR + x M). Assume that the c^ 1 and c have a common lower bound 

//el. 

Then, for all a, (3 G M + satisfying < a < (3 < oo, S(vq ,c^) converges to S(vq,c) in 
L 2 ([a, 0] xl). 

Proof . Let us introduce the intermediate S(t;o,c' n '). Then by the preceding Proposition 

|i% Q ,c^) - 5(^ n) , C W)|| L2(K/3]xR) < |bo - 4 n) |Us (R ) f / e~ 2 ^ ^ 
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and this term converges to 0. 

The second term S(vq, c^) — S(v , c) is more delicate. Let 



I n {t,x) 



Vol 



2tw{\) + x ) e -*Jo(*+VS«i,to) a d. 



c w 



e -t c^(t{l~s),V2tw{s)+x)ds _ e -tfic(t(l-s),V2iw(s)+x)ds^ ^ 

so that we may write 

\\S{v Q , c (n) ) - S(v , c)\\ 2 L 2 MxR) = [ [ I n (t,x) 2 dt dx. 

J a JR 

We shall prove that the sequence (I n (t, x)) n converges to when n goes to infinity and that 
it is uniformly bounded by a function g(t, x) belonging to L 2 ([a, 0] xl). Then, thanks to 
the dominated convergence Theorem, | \S(v , c^ 71 ') — S(v , c)\\\a(\ a m x jg\ will converge to 0. 

As usual, the first exponential term of I n (t,x) is smaller than 1. Both —t c^ n \t(l — 
s), y/2tw(s) + x)ds and —t f Q c(i(l — s), y/2tw(s) + x)ds being smaller than — tfj,, 

-t Jg cW(t(l-s),V2iw{s)+x)ds _ -t Jg c( y t(l-s),V2tw(s)+x^ds 



< 



t! (c< n) -c)(t(l-s), V2tw(s) + x)ds 
Jo 



-tp, 



and we can estimate \I n \ by 



I n (t,x)\ < te-* [ \v (V2tw(l) + x)\ f (c^ -c)(t(l-s),V2tw(s) + x) 

JC'w Jo 



dsdmw 



M n {t,x) 

Thanks to (J4]), M n (t,x) can be written as an integral on 1 



M n (t,x) 



[ f \vo(V2k 2 + x)\ {c^-c){t(l-s),V2t^+x) f s ,i(Zi,&)d£id&ds 
Jo Jr 2 



where / Sj i is the gaussian density. The change of variables 

u = y/2t£,2 + x, v = t(l — s), w = V2t £i + x 

gives 



1 



Mn(t,x) =— . 

47Tt J R 3 y/ V (t - V) 



ltn t ](v) , / \ . (w-x) 2 (u-w) 2 

'(c w - c)(v, w) e 4 (*-") e ^ \v Q (u)\ dudvdw. 



As 



(if — in) 



e *» \v (u)\ du < \\v \\ L 2 m (27rv) 



1/4 
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M n (t,x) is smaller than 



1 II II /o \l/4n (Will f 1 [0,*]( U ) 1/9 A (W-X) 2 



dvdw. 



The integral appearing in the square root converges and does not depend on n, which 
shows that M n (t,x) and I n (t,x) go to 0. 
Now for the uniform estimate. Clearly 



|J„(*,x)|< / \v (V2tw(l) +x)\(e- tlM + e-^)dm w := g(t,x). 



The function g is in L 2 ([a, (3} x M.) since 

/ f g(t,x) 2 dxdt < [ 4e- 2 ^ /Y / 

Ja JR \Jc 

< [ 4e~ 2tfl [ [ \v (y/2tw(l) + x)\ 2 dm w dxdt 

J a JR JCw 

< f\e- 2t » [ [ \v (V2tw(l) + x)\ 2 dx dm w dt . 

J a JCw JR 



p r rP r / r \ 2 

v Q (y2tw(l) + x)\dmw I dxdt 



The change of variables £ = y/2tw(l) + x allows to write 

/ / g(t,x) 2 dxdt < \\v \\ 2 L 2( R) < / Ae~ 2t>1 dt < oo, 

Ja JR Ja 

which concludes the proof. □ 

Any vq G L 2 (M) can be approximated (in L 2 (IR)) by a sequence (vq) of functions 
satisfying the hypotheses of Theorem [5j Similarly, any function c belonging to L 2 (]0, oo) x 
R) n L°°(]0,oo) x E) is the limit (in L 2 (]0,oo) x R)) of a sequence (c (n) ) n of functions 
satisfying the hypotheses of Theorem [5] Moreover we can suppose this sequence to be 
bounded in L°°(]0, oo) x M). Then, according to Theorem El v n := S(vq\ c^) is a solution 
of 

^(t,x)-^(t,x) + (x 2 + c^(t,x))v n (t,x) = on ]0,oo[xM 
at ox z 

v n (0,x) = v^\x) . 

Let tp be a smooth test function on ]0, oo[xR. Integrations by part give 

dcp d 2 ip 2 (n) 

< v n , — > - < v n , — T > + < x v n , <p > + < & >v n , <p >= 0, 



dt ' dx 
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No- 



where the brackets stand for L 2 (]0,oo[xR) products. As v n converges to S(vq,c) in any 
L 2 ([a,f3] x R), we obtain 

dip 2 
< S(v , c), — > - < S(v , c), — > + < x 2 S(v Q , c), ip > + < cS(v , c), if >= 0. 

Moreover, v n (0, •) = v^ n ' converges to v , which gives the equality S(v , c)(0, •) 
Hence S(vq, c) is a solution of ([1]) in a weak sense. 

Appendix A. Computation of some Wiener integrals 
Another consequence of Theorem [1] is the following Proposition : 
Proposition 12. Suppose Vq satisfies, for all positive t and real x 



[ (2tt sh(2t))- 1/2 exp (- 



lch(2t) 2 lsh(2t) 2 



2sh(2ff 2ch(2t) 
Then for allt>0 and x G R we can write 



vo[y + 



X 



ch(2t) 



dy < oo. 



J Vo(V2t w(l) + x) exp ^— t J (x + y/2t w(s)j ds^j dmwiw) 

I (2vr sh(2t))~ 1/2 exp 
Jr 



lch(2t) 2 lsh(2t) 2 
2sh(2t) y ~ 2ck(2t) X 



vo\y + 



X 



ch(2t) 



dy 



Proof. Suppose vq belongs to L 2 (R). Then Theorem CD shows that, as the perturbation 
c is equal to 0, the left hand side is the solution of the heat equation for the harmonic 
oscillator, with initial condition vq. The right hand side is the solution UfVQ of the same 
problem, given by Mehler's formula. The equality follows. 

When vq is not in L 2 (R), the equality holds for v^ifn where <f n is a convenient truncature. 
Then the assumption on vq allows to use the Theorems of dominated and of monotone 
convergence. □ 

When Vq(x) = l,x or x 2 , it is easy to compute U t vo and to deduce the following 
equalities from these computations : 



Corollary 13. For vq = 1 we obtain 



k(t,x) :- 
The case vq(x) 



1 



e -tti( X+ V2i «,(.)) ds d ( w ) = 

c w v^M2t) 

x gives 



w(l) exp —t 



Cw 



X 



2t wis 



ds ) dmwiw 



1 sh(2t) 
6XP ^~2ch(2t) 

x(l - ch(2t)) 



x 



and for Vq(x) = x 2 we get 



/ u>(l) 2 exp(— t / (x+V~2t w(s)) 2 ds) dmw\ 
Jc w Jo 



w) — — 
' 2t 



,— k(t,x) 
f2i ch(2t) 

1 - ch(2t))V sh(2t)\ 



ch 2 (2t) 



+ 



ch(2t) J 



k(t, x) . 
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As the perturbation c is equal to 0, [3] and [6] prove that, under certain conditions, the 
l.h.s. is the solution of the heat equation problem for H. Nevertheless these integrals are 
not mentioned explicitly in the literature, to the author's knowledge. 
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